Electronic Raman response in anisotropic metals 
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Using a generalized response theory we derive the electronic Raman response function for metals 
with anisotropic relaxation rates. The calculations account for the long-range Coulomb interaction 
and treat the collision operator within a charge conserving relaxation time approximation. We 
extend earlier treatments to finite wavenumbers (|q| <C acf) and incorporate inelastic electron- 
electron scattering besides elastic impurity scattering. Moreover we generalize the Lindhard density 
response function to the Raman case. Numerical results for the quasiparticle scattering rate and 
the Raman response function for cuprate superconductors are presented. 



Introduction 

After the discovery of high-T c superconductors, some- 
thing of a stir has been caused by the possibility of inves- 
tigating these systems by Raman scattering experiments. 
Through its polarization dependence this spectroscopic 
method allows for detecting anisotropics in quantities 
like the superconducting gap and the scattering rates 
characterizing the transport properties. While the gap 
anisotropics in high-T c cuprates, as seen in Raman ex- 
periments, are theoretically well understood in terms of 
d y 2_ y 2 pairing 1 , the important quasiparticle collision ef- 
fects are much less studied. Early attempts included con- 
stant scattering rates 2 , elastic scattering processes 3,4 or 
described inelastic scattering processes within the Nearly 
Antifcrromagnctic Fermi Liquid (NAFL) model 5 . 

In this paper we investigate the general situation of 
a normal metal by developing a theory of the electronic 
Raman effect at finite wavenumbers q =/= within the 
RPA response theory and generalize it to include colli- 
sion effects. A superposition of elastic and inelastic scat- 
tering processes is considered with k-dependent relax- 
ation rates and deviations from Matthicssen's rule are 
studied. Having in mind quasi 2D-systcms like high-T c 
cuprate superconductors, we present a numerical analysis 
of the inelastic scattering rates and the Raman response 
functions employing the FLuctuation EXchange (FLEX) 
approximation which treats the spin fluctuation limited 
transport and Cooper-pairing on the same level. 
Transport theory and generalized Raman 
response: normal state 

We consider a normal metal in which the electronic states 
are characterized by a momentum Kk, an energy disper- 
sion ek = £k + fi, (we set the lattice constant to unity) 

6k = t [—2 cos(k x ) — 2 cos(ky) + 4Bcos(k x ) cos(fc 3/ )] , (1) 

with [i being the Fermi energy, the group velocity Vk = 
(l/7i)Vkek, an inverse effective mass tensor M-~ 1 (k) = 
d 2 e\ i /h 2 dkidkj 7 and an equilibrium Fermi distribution 
nk with derivative ip-^ = —dn^/d^. Next we focus on 
an external perturbation appropriate for a treatment of 
the electronic Raman response within the effective mass 
approximation: 6-8 



C/ k cxt = me s ■ M-!(k) • e 1 • -i-IA^AI . (2) 




Here, e 1 and e s are the unit vectors of the incident and 
scattered light, respectively, and A denotes the vector 
potential. Then, the response of the electronic system to 
this perturbation is described within the the quasiclassi- 
cal limit of the kinetic equation: 

u5n k - q • v k h k = i ^ 5I£ , (3) 

v— e,i 

where we use the definition 

h k = 8n k + <p k [U£ xt + V(q)dn 1 ] (4) 

and V(q) = 47re 2 /q 2 is the Fourier transform of the long 
range Coulomb interaction. Physical observables are the 
generalized response functions 

8n a = ^ a p Sn p , (5) 

where a p is the vertex which describes the coupling of 
5n a to the external perturbation potential. The collision 
integrals for elastic (y =e) and inelastic [y =i) scattering 
have the form (conserving relaxation time approxima- 
tion) 

per 

gk P ~yk^A 6 . (6) 

b l^pa ^PV p 

The scattering parameters allow a classification of the 
macroscopic moments Srib (defined through Eq. (5)) into 
conserved (A^ = 1) and nonconserved (A^ < 1) quanti- 
ties. In what follows we will, for the sake of simplicity, 
restrict ourselves to the case of charge conservation, i.e. 
\\ = 1, XI = V6 ^ 1. While the elastic scattering 
rate T k is constant, the inelastic part of the scattering 
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rate has a strong frequency dependence that reflects 
characteristic lifetime effects, which we discuss below. 

In Fig. 1 we show the inelastic scattering rate 
r k (w) for hole-doped high-T c cuprates at optimal dop- 
ing (x = 0.15) for various temperatures using the FLEX 
approximation 9-12 . The upper curves correspond to In- 
directions near (tt,0) ('hot spots'), the lower ones refer 
to k near the diagonal ('cold spots'). For u> — > we 
find that the inelastic scattering rate at the hot spot is 
almost three times larger than at the cold spot (solid 
line). Physically speaking, the scattering rates can 
be understood in terms of scattering of quasiparticles 
on spin fluctuations (paramagnons) which are enhanced 
near the hot spots. In the normal state we also find at 
small frequencies that the scattering rates decrease with 
decreasing temperature. For T < T c a rearrangement 
of spectral weight occurs reflecting the w-dependence of 
the superconducting d a .2_. y 2-wave gap A(uj) which is cal- 
culated self-consistently 13 . As expected, the maximum 
is seen at oj ~ 3A /?i. In the high-frequency limit, T l k 
varies linearly with to for all k-dircctions. This is in ac- 
cordance with both the Marginal (MFL) 14 and Nested 
(NFL) 15 Fermi liquid picture. 



> 

s 



3, 

— M 

u 



160 



100 
80 
60 
40 
20 




x = 0.15 




'hot spot' 






'cold spot' 



20 40 60 80 

co (meV) 



1 00 



120 



FIG. 1. Calculated inelastic scattering rate r k versus uj for 
the one-band Hubbard model using U = 4£ at optimal doping 
(x = 0.15) for temperatures T = 2T C (solid lines), T = 1.05T C 
(dashed lines), and T = 0.5T C (dashed-dotted lines). 

Let us turn to the overdoped case (x = 0.22) where we 
focus on the normal state. Thus, we show in Fig. 2 for 
the same temperatures and directions as in Fig. 1. Most 
importantly, the scattering rates become less anisotropic 
at small frequencies. This is in agreement with Raman 
scattering experiments where the static scattering rate is 
extracted for different high-T c cuprates for different po- 
larisations as a function of the doping concentration 16 ' 17 . 
Generally speaking, the anisotropy of the scattering rate 
reflects the anisotropy of the spin fluctuations. Those 
become less pronounced and more isotropic in the over- 
doped regime. We still find a reduction of T l k with de- 
creasing temperature and the linear high-frequency be- 
havior remains. 



Next, we turn to the generalized theory of Raman scat- 
tering. The full response function C a b for general vertices 
a and b reads within RPA (see (Eq. 2)) 
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FIG. 2. Inelastic scattering rate F k versus u> for the over- 
doped case (x = 0.22) using the same notation as in Fig. 1. 
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/ denotes the corresponding irreducible interaction in a 
symbolic notation. Note that the decomposition into 
transverse and longitudinal parts is a general structural 
property of the respose formalism. For the special case 
of Raman scattering (a, b = 7k) one has to consider the 
Lindhard response and (1 — ILu) is the dielectric func- 
tion e. Then, the full Raman response is of the form 
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Here, M* b is the generalization of the Lindhard function 
due to the inclusion of collision effects: 
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The quantities are generalized diffusion tensors 
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S* b are the collision-limited Raman response functions 
which have a finite q — > limit: 
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and the objects (* b describe the mixing of clastic and 
inelastic scattering processes and hence deviations from 
Matthicssen's rule within the conserving relaxation time 
approximation: 
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Finally, e = e(q, u) = 1 — V(q)Mj* 1 (q, cj) is the dielectric 
function of the electronic system mentioned earlier. It 
has thus been shown that the mixing terms oc £* b in Eq. 
(8) of the two separate scattering mechanisms occur in 
the Lindhard function M* h (q, us) = 0(q 2 ) and in terms 
which are screened oc e^ 1 by the long range Coulomb in- 
teraction. For practical applications of the result shown 
in Eq. (8) to the cuprate systems, say, where q — > and 
e rj 10 4 may be assumed, one may use the q — > 0, e — > oo 
limit of L 77 and the contributions from mixing oc £* h are 
hence irrelevant. Therefore, the scattering mechanisms 
can be linearly combined, leading to the exclusive occur- 
rence of r* = r p + T p in the Raman response function. 
Note that the Raman response function Im £ 77 
limit q — ► reduces to Eq. (12) of Ref. 3 where (r£ 
corresponds to an average of T*. Hence, this scattering 
rate can be interpreted as the width of the lineshape of 
Fig. 3 in Ref. 3. Next we focus on the numerical analysis 
of Im £ 77 for a spin-fluctuation based model particularly 
for the inelastic scattering rate T^. 
Numerical calculation of the Raman response 
Using normal and anomalous Green's functions, G and 
F, we calculate the full Raman response function via 

/oo 
duj' [/(a/) - f(u' + u)] 
-OO 

x 7( k ' w ) [ N (^ J + u)N(k, a/) 

k 

- A(k,Lu' +cu)A(k,Lu')] 1 (k) , (13) 

where, N — G" and A = F" are the spectral functions 
(i.e. the k- resolved density of states) of the corresponding 
quasiparticles calculated within the FLEX approxima- 
tion for the one-band Hubbard model 9 ' 18 . The quantity 
7(k, u>',w) denotes the renormalized Raman vertex that 



will be specified below. The bare Raman vertices for the 
different polarization symmetries considered here, Bi g , 
£?2g, read within the effective mass approximation: 

7Bi s = t [cos(fc x ) - cos(fc y )] ,7s 2g = 4tBsm(k x ) sm(k y ). 

Here, t is the nearest neighbor and t' = —Bt (with 
B = 0.45) is the next-nearest neighbor hopping energy of 
the tight-binding band introduced in Eq. (1). Thus Ra- 
man scattering in B\ g symmetry mainly probes the 'hot 
regions' while B^ g symmetry probes the 'cold regions' of 
the Brillouin zone. 

It has been shown earlier that the FLEX approxima- 
tion yields a phase diagram T{x) for cuprates which is in 
fair agreement with experiment. One finds a d K 2_ y 2-wave 
superconducting order parameter and all characteristic 
temperature scales 19 . Recently it has also been demon- 
strated that the resonance peak below T c in neutron scat- 
tering data and angle-resolved photoemission (ARPES) 
experiments which measure the spectral density entering 
in Eq.(13) can be well described 20,13 . 

While no vertex corrections are considered using the 
kinetic equation approach, they become important if 
the FLEX approximation is employed. The main con- 
sequence is that the Raman vertex becomes frequency 
dependent reflecting the retardation effects related to 
the Cooper-pairing mechanism. In order to calculate 
the vertex function 7„ we employ the Nambu notation 
(//,i/ = 0,l,2,3) 21 
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with q = q, iv m ; k = k, iuj n ; J^k = T Sk> and sat " 

isfy the Ward identity £^ <? M 7 M (fc + q,k) = r 3 G^ 1 (k) - 
G^ 1 (k + q)r 3 . The ladder approximation for the vertex 
function yields 

%(k + q,k) =7 M (fc + (?,fc) 

+ I T -oG(fc / + q)%{k' + q, k')G(k')T Q P s (k - k') 

k' 

+ T 3 G(k' + q)%(k' + q, k')G(k')T 3 P c (k - k')} . (15) 

P s and P c refer to the spectral function of the spin and 
charge susceptibility within RPA, respectively, that are 
defined in Ref. 9. Finally, the full Raman vertex 7 satis- 
fies the following integral equation: 

j(k + Q,k)= 1 (k) 



j(k + q + Q, k + q) 



(16) 



and can be derived from Eq. (14) by replacing 70 and 70 
by 7T3 and 7T3. 

Results for £ 77 and discussion 

In Fig. 3 we show our results for the full Raman re- 
sponse function Eq. (8), £ 77 , for the B\ g polarization, 
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which is proportional to the measured Raman intensity, 
at optimal doping. Both a linear increase and a contin- 
uum are found in the limit of low and high frequencies, 
respectively. This is a direct consequence of the inelastic 
scattering rate discussed in Fig. 1. In particular, we find 
an increasing initial slope with decreasing temperature. 
Below T c we obtain the usual pair-breaking feature 
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FIG. 3. Full Raman response function £ 77 (q —> 0,lo) ver- 
sus lo at optimal doping (x = 0.15) for B\ g polarization for 
the same temperatures as in Fig. 1. 

accompanied by a suppression of spectral weight for 
u> < 2Ao/7i. For low frequencies we also find a power 
law £ 77 oc (oj/2Ao) 3 that is charcteristic for a d x i_ y i- 
wave gap in the clean limit. Note that the pair breaking 
peak is finite and renormalized by inelastic quasiparticle 
scattering processes. 

Figure 4 shows the Raman intensity for i?2 9 
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FIG. 4. Full Raman response function £ 77 (q —> 0,lj) ver- 
sus uj at optimal doping (x — 0.15) for B-2 g polarization for 
the same temperatures as in Fig. 1. 

polarization in the normal state. We again find the 
high-frequency continuum and an increasing slope of the 
Raman response function with decreasing temperature. 
This agrees with the NAFL picture if vertex corrections 
are considered 5 . 



Summary 

In summary we have reconsidered the electronic Raman 
response with special emphasis on the role of quasiparti- 
cle scattering processes. In the normal state a superpo- 
sition of elastic and inelastic scattering events described 
by k-dependent scattering rates leads to a deviation from 
Matthicssen's rule for the corresponding Raman response 
functions. These deviations, however, are found to be 
eliminated by the presence of the long-range Coulomb 
interaction. We have furthermore studied in detail, the 
anisotropy of the inelastic scattering rate at opti- 
mal doping, and its decrease at higher doping levels to- 
wards the overdoped case. The resultant electronic Ra- 
man response functions that include the effects of inelas- 
tic scattering involving spin fluctuations reproduce the 
behavior of the normal state characteristic of the phe- 
nomenological MFL picture that has been studied in ear- 
lier publications 22 . In the superconducting state the in- 
elastic scattering rate becomes gapped due to a frequency 
dependent d x i _j,2-wave gap, calculated self-consistcntly, 
yielding the characteristic shape of the Raman response. 1 
Acknowledgements 

D.M. thanks the Walthcr-Mcissncr-Institute for hospital- 
ity and INTAS (project 01-0654) for financial support. 



1 T.P. Devereaux and D. Einzel, Phys. Rev. 51, 16336 
(1995); Phys. Rev. 54, 15547 (1996). 

2 L.A. Falkovsky and S. Klama, Physica C 172, 242 (1990). 

3 A. Zawadowski und M. Cardona, Phys. Rev. 42, 10732 
(1990). 

4 D. Einzel and C. Schuster, Czechoslovak J. of Physics 46, 
993 (1996). 

5 T. P. Devereaux and A. P. Kampf, Phys. Rev. 61, 1490 
(2000); Phys. Rev. 59, 6411 (1999). 

6 P. A. Wolff, Phys. Rev. 171, 436 (1968). 

7 M. C. Krantz and M. Cardona, J. Low Temp. Phys. 99, 
205 (1995); T. Strohm and M. Cardona, Phys. Rev. B 55, 
12725 (1997). 

8 D. Manske, C. T. Rieck, R. Das Sharma, A. Bock, and D. 
Fay, Phys. Rev. B 56, R2940 (1997). 

9 See, for example, D. Manske, I. Eremin, and K. H. Benne- 
mann, Phys. Rev. B 67, 134520 (2003). 

10 N. E. Bickers, D. J. Scalapino, and S. R. White, Phys. Rev. 
Lett. 62, 961 (1989); N. E. Bickers and D. J. Scalapino, 
Ann. Phys. (N.Y.) 193, 206 (1989). 

11 T. Dahm and L. Tewordt, Phys. Rev. Lett. 74, 793 (1995). 

12 M. Langer, J. Schmalian, S. Grabowski, and K.H. Benne- 
mann, Phys. Rev. Lett. 75, 4508 (1995). 

13 D. Manske, I. Eremin, and K. H. Bennemann, Phys. Rev. 
Lett. 87, 177005 (2001). 

14 C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abra- 
hams, and A. E. Ruckenstein, Phys. Rev. Lett. 63, 1996 
(1989). 

15 J. Ruvalds, CT. Rieck, S. Tewari, J. Thoma, and A. Vi- 
rosztek, Phys. Rev. B 51, 3797 (1995). 

16 For a systematic study, see M. Opel, R. Nemetschek, C. 
Hoffmann, R. Philipp, P. F. Miiller, R. Hackl, I. Tutto, A. 



4 



Erb, B. Revaz, E. Walker, H. Berger, and L. Forro, Phys. 
Rev. B 61, 9752 (2000). 

F. Venturini, M. Opel, T. P. Devereaux, J. K. Freericks,I. 
Tutto, B. Revaz, E. Walker, H. Berger, L. Forro, and R. 
Hackl, Phys. Rev. Lett. 89, 107003 (2002). 
Note that the NAFL approach of Ref. 5 is restricted to the 
normal state and thus no spectral function A occurs. 
D. Manske, I. Eremin, and K. H. Bennemann, Phys. Rev. 
B 62, 13922 (2000); D. Manske, T. Dahm, and K. H. Ben- 
nemann, Phys. Rev. B 64, 144520 (2001). 
D. Manske, I. Eremin, and K. H. Bennemann, Phys. Rev. 
B 63, 054517 (2001). 
Y. Nambu, Phys. Rev. 117, 648 (1960). 
T. Staufer, R. Hackl, and P. Miiller, Solid State Comm. 
75, 975 (1990). 



5 



